Skip to main content

4.2 THRUST CHAMBER PERFORMANCE PARAMETERS

The significance of the parameters which express or influence the efficiency of thrust chamber operation has been discussed in section 1.3 of chapter I. Before discussing the details of actual thrust chamber design, the following summarizes these parameters and illustrates their use for design calculations by applying them to the engine systems of the Alpha vehicle which was discussed in chapter III.

Specific Impulse, IS(sec)I_{S}(\mathrm{sec})

From equations 1-31 and 1-31c:

(Is)tc=FW˙tc=cCfg\left(I_{\mathrm{s}}\right)_{\mathrm{tc}}=\frac{F}{\dot{W}_{\mathrm{tc}}}=\frac{c^{*} C_{f}}{g}

The specific impulse figure indicates the overall quality of the thrust-chamber design. As was learned earlier, it shows how much thrust is generated for what "running propellant expenditures."

Characteristic Velocity, c* (ft/sec)

From equation 1-32a:

c=f(γ,R,(Tc)ns)\begin{equation*} c^{*}=f\left(\gamma, R,\left(T_{c}\right)_{\mathrm{ns}}\right) \tag{4-1} \end{equation*}

Assuming that the propellant and mixture ratio selection has been made, it can be expected that the gas properties ( γ,R\gamma, R ) will fall into a known band. From there on, cc^{*} almost entirely depends on the temperature of the gases. Obviously, this temperature has a theoretical maximum for a selected propellant combination. How close to this maximum the chamber will operate depends on influences discussed for mixture ratio in chapter II. Figures 43,44,454-3,4-4,4-5, and 464-6 illustrate this. It is seen that cc^{*} peaks at combustion temperatures somewhat lower than maximum. Other considerations, such as bulk densities, which affect vehicle tank sizes, may cause further adjustments of the mixture ratio for optimum overall vehicle performance. Within these boundaries, the quality of the combustion process greatly depends on the design efficiency of a thrust chamber assembly, in particular the injector.

Thrust Coefficient, CtC_{t} (dimensionless)

From equation 1-33a:

Cf=f(γ,ϵ,pa)\begin{equation*} C_{f}=f\left(\gamma, \epsilon, \mathrm{p}_{\mathrm{a}}\right) \tag{4-2} \end{equation*}

Let us assume that the performance of energy generation through the combustion process, the effects of which were just summarized for cc^{*}. has been determined. Then, with a given gas property ( γ\gamma ), the performance of the remaining thrust generating functions of the thrust chamber, essentially those of the divergent nozzle, will depend on the nozzle geometry (mainly ϵ\epsilon which determines the pressure ratio (pe/(pc)ns)\left(p_{e} /\left(p_{c}\right)_{n s}\right) and the ambient pressure ( pap_{\mathrm{a}} ).

Performance Calculation

In actual design practice, the calculation of the thrust chamber performance is based on theoretical propellant combustion data and the application of certain correction factors as explained in chapter I. The theoretical propellant combustion data are derived from thermochemical computations which equate the heat of reaction of the propellant combination and the rise in enthalpy of the combustion gases. Typical propellant combustion data at frozen composition are presented in figures 434-3 through 4-6. For given propellant combinations and chamber-nozzle stagnation pressures (Pc)ns\left(P_{c}\right)_{\mathrm{ns}}, the values for the combustion gas temperature (Tc)ns\left(T_{c}\right)_{\mathrm{ns}}, molecular weight M1M_{1}, and specific heat ratio γ\gamma are plotted against the O/F\mathrm{O} / \mathrm{F} mixture ratio rwr_{w}. Performance correction factors are determined by the theoretical assumptions and from earlier test data, as well as from the chosen design configurations. Typical performance calculation methods have

Figure 4-3.-Theoretical O2/RP1\mathrm{O}_{2} / R P-1 combustion data (frozen composition), (pc)ns=1000psia\left(p_{c}\right)_{\mathrm{ns}}=1000 \mathrm{psia}.

Figure 4-4.-Theoretical O2/H2\mathrm{O}_{2} / \mathrm{H}_{2} combustion data (frozen composition), (pc)ns=800psia\left(p_{c}\right)_{\mathrm{ns}}=800 \mathrm{psia}.

Figure 4-5.-Theoretical F2/H2F_{2} / H_{2} combustion data (frozen composition), (pc)ns=100psia\left(p_{c}\right)_{\mathrm{ns}}=100 \mathrm{psia}.

Figure 4-6.-Theoretical N2O4/N2H4\mathrm{N}_{2} \mathrm{O}_{4} / \mathrm{N}_{2} \mathrm{H}_{4} combustion data (frozen composition), (pc)ns=100psia\left(p_{c}\right)_{\mathrm{ns}}=100 \mathrm{psia}.

been demonstrated earlier by sample calculation (1-3). The following sample calculations illustrate the more specific approaches.

Sample Calculation (4-1)

Determine the design values of c,Ctc^{*}, C_{t}, and (Is)tc\left(I_{\mathrm{s}}\right)_{\mathrm{tc}} for the engine thrust chambers of the stages of the hypothetical Alpha vehicle, with the following assumed design parameters: (a) 750 K A-1 Stage Engine:

Propellants, LO2/RP1\mathrm{LO}_{2} / \mathrm{RP}-1; thrust chamber O/F mixture ratio, 2.35; (Pc)ns ,1000\left(P_{c}\right)_{\text {ns }}, 1000 psia; propellant combustion data, figure 4-3; nozzle expansion area ratio, ϵ=14\epsilon=14. (b) 150K150 K A-2 Stage Engine:

Propellants, LO2/LH2\mathrm{LO}_{2} / \mathrm{LH}_{2}; thrust chamber O/F mixture ratio, 5.22;(pc)ns ,8005.22 ;\left(p_{c}\right)_{\text {ns }}, 800 psia; propellant combustion data, figure 4-4; nozzle expansion area ratio, ϵ=40\epsilon=40.

Solutions

(a) A-1 Stage Engine:

From figure 4-3 for LO2/RP1\mathrm{LO}_{2} / \mathrm{RP}-1 at (Pc)ns=1000\left(P_{c}\right)_{\mathrm{ns}}=1000 psia and a mixture ratio of 2.35 , the following values are derived for the chamber product gases:

(Tc)ns=6000F or 6460R,π=22.5lb/mol,y=1.222\begin{aligned} & \left(T_{c}\right)_{\mathrm{ns}}=6000^{\circ} \mathrm{F} \text { or } 6460^{\circ} \mathrm{R}, \\ & \quad \pi=22.5 \mathrm{lb} / \mathrm{mol}, \quad y=1.222 \end{aligned}

Substitute into equation (1-32a): Theoretical c=32.2×1.222×6460×1544/22.50.7215c^{*}=\frac{\sqrt{32.2 \times 1.222 \times 6460 \times 1544 / 22.5}}{0.7215}

=5810ft/sec=5810 \mathrm{ft} / \mathrm{sec}

This value for cc^{*} can also be derived from figure 4-3.

For a good combustion chamber and injector design, the cc^{*} correction factor for LO2/RP1\mathrm{LO}_{2} / \mathrm{RP}-1 and frozen composition will be about 0.975 .

Design c* =5810×0.975=5660ft/sec=5810 \times 0.975=5660 \mathrm{ft} / \mathrm{sec} For y=1.222,ϵ=14y=1.222, \epsilon=14, a theoretical vacuum CfC_{f} value of 1.768 can be derived from figure 1-11:

 Theoretical Cf at sea level =(Cf)vacϵPa(Pc)ns=1.76814×14.71000=1.562\begin{aligned} \text { Theoretical } C_{f} \text { at sea level } & =\left(C_{f}\right)_{\mathrm{vac}}-\frac{\epsilon P_{a}}{\left(P_{c}\right)_{\mathrm{ns}}} \\ & =1.768-\frac{14 \times 14.7}{1000}=1.562 \end{aligned}

Sea level CfC_{f} can also be calculated using equation (1-33a), with the aid of equation (1-20).

With effective nozzle contour design, an overall CfC_{f} correction factor of 0.98 for LO2/RP1\mathrm{LO}_{2} / \mathrm{RP}-1 frozen composition can be used.

Design sea level Cf=1.562×0.98=1.531C_{f}=1.562 \times 0.98=1.531 From equation (1-31c): Design sea level (Is)tc=5660×1.53132.2=270sec\left(I_{s}\right)_{\mathrm{tc}}=\frac{5660 \times 1.531}{32.2}=270 \mathrm{sec} (b) A-2 Stage Engine:

From figure 444-4 for LO2/LH2\mathrm{LO}_{2} / \mathrm{LH}_{2} at (Pc)ns=800\left(P_{c}\right)_{\mathrm{ns}}=800 psia and a O/F\mathrm{O} / \mathrm{F} mixture ratio of 5.22 , the following values are derived for the chamber product gases: (Tc)ns=5580F\left(T_{c}\right)_{\mathrm{ns}}=5580^{\circ} \mathrm{F} or 6040R6040^{\circ} \mathrm{R},

π=12lb/mole,y=1.213\pi=12 \mathrm{lb} / \mathrm{mole}, \quad y=1.213

Substitute into equation (1-32a): Theoretical c=32.2×1.213×6040×1544/120.717=7670ft/secc^{*}=\frac{\sqrt{32.2 \times 1.213 \times 6040 \times 1544 / 12}}{0.717} =7670 \mathrm{ft} / \mathrm{sec}

Based on experimental data, a cc^{*} correction factor can be assumed for the LO2/LH2\mathrm{LO}_{2} / \mathrm{LH}_{2} frozencomposition data of about 0.975 .

 Design c=7670×0.975=7480ft/sec\text { Design } c^{*}=7670 \times 0.975=7480 \mathrm{ft} / \mathrm{sec}

For y=1.213,ϵ=40y=1.213, \epsilon=40, a theoretical vacuum CfC_{f} value of 1.876 can be derived from figure 1-11. CfC_{f} can also be calculated using equations (1-33a) and (1-20). With effective nozzle contour design, an overall CfC_{f} correction factor value of 1.01 can be used for LO2/LH2\mathrm{LO}_{2} / \mathrm{LH}_{2} frozen-composition data:

Design vacuum Cf=1.876×1.01=1.895C_{f}=1.876 \times 1.01=1.895 From equation (1-31c): Design vacuum (Is)tc=7480×1.89532.2=440sec\left(I_{\mathrm{s}}\right)_{\mathrm{tc}}=\frac{7480 \times 1.895}{32.2}=440 \mathrm{sec} The reader should perform his own calculations for the A-3 and the A-4 engines, with the aid of tables 3-4 and 3-5, and figures 4-5 and 4-6.